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A preliminary version of the compressible second-order closure model that will be con- 
sidered for use in the National Aero-Space Plane Project (NASP) is presented on the pages 
to follow. The proposed model requires the solution of modeled transport equations for 
the Favre-averaged Reynolds stress tensor and the turbulent dissipation rate in addition to 
the corresponding mean continuity, momentum and energy equations. In this preliminary 
model, the Reynolds heat flux and the mass flux terms are modeled by a simple gradient 
transport hypothesis (the same is true for the turbulent diffusion terms in the Reynolds 
stress transport equation). A model recently developed by Speziale, Sarkar and Gatski 
(1989) for the pressure-strain correlation (which more properly accounts for rotational 
strains) is used for the deviatoric part of the pressure gradient-velocity correlation. This 
early version of the model neglects correlations involving fluctuating dilatational terms. 
The modeled dissipation rate transport equation has been formulated to be consistent 
with Rapid Distortion Theory for a compressed isotropic turbulence. In order to facilitate 
the early implementation of the model into existing computer codes, Van Driest damping is 
proposed for integration to the wall. Ultimately, this will be supplanted by asymptotically 
correct low Reynolds number corrections to the models in order to obtain a more accurate 
calculation of wall flow properties. It is envisioned that as this modeling effort progresses, 
transport equations for the Reynolds heat flux, the mass flux, and possibly the rms density 
and temperature fluctuations may be added. Likewise, the modeling of dilatational effects 
(e.g., a model for the pressure-dilatation correlation) and the possible addition of an alter- 
nate specification for the turbulence length scale are being investigated with J. L. Lumley 
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(Cornell University) and T. B. Gatski (NASA Langley Research Center). A systematic 
program for the testing of this model which makes use of the results of both physical and 
numerical experiments is currently being developed in collaboration with T. B. Gatski and 
A. Kumar (NASA Langley Research Center). 
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Preliminary Compressible Second-Order Closure Model 


Nomenclature 


p = mass density 
Vi = velocity vector 
p = pressure 
T = temperature 

C v = specific heat (constant volume) 

Qi = Favre-averaged Reynolds heat flux 

Tij = Favre-averaged kinematic Reynolds stress tensor 

k = turbulent kinetic energy 

a ij = viscous stress tensor 

p = dynamic viscosity 

J? = ideal gas constant 

$ = viscous dissipation rate 

e = turbulent dissipation rate 

k = thermal conductivity 

a,( ) = a( )/at 

( ),*' — )/^ x * 


For any flow variable / : 

/ = / + /" , (J=\im±J Q f{x,t)dt) 

f = ~f + f' , (/= ^ = Favre average ) 

Mean Continuity equation: 


d t (p) + (/ft*),* = 0 


( 1 ) 
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Mean Momentum equation: 


dt(pVi) + (pViV^j = - Pi + - {pTij) tj (2) 

where 

2 

Vij = --Wkjkhi + + V 3,i) 

2__ f „ . 

- - ^Vk,kOij + mK j + v it i) 

Tij = v'u'- , p = pRT 

Mean Energy equation: 

dt(pC v T) 4- {pvkC v T) t k = —pvkji + $ + («rjt),k — Qk,t (3) 

where 

$ = + a ij V lj 

- VijVij + pe 

WkJ - PVk.k 
= 7c T k 

Q k = JfC v vP> ~ -JSC^—Tk 

e&T 

The model parameters required in the above expressions are 

Oj. = 0.7 


and 


= 0.09[l — exp(— y + /25)] 
y + = pu T y/ji 
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Favre-averaged Reynolds stress transport equation: 


The exact transport equation for the Favre-averaged Reynolds stress r,y is as follows 


d t {pTij) + (vkPTij),k = -pTikVj,k - pTjkVi,k - Cijk'k - II, y - €,y + 

-%Pj - + tyik,k + (vy[ k + v'jo'ik) ,k (4) 


where 


Cijk 

n„ 


= P V 'i V j V k + ~P"< 6 i} 


m = 


*<fi + 

°>",k + °%<,k 


The models that are proposed in order to close (4) are as follows 


2 _ 

£ 


V WA - -ff. T (rnjk + nkj + r iV ) 

pX = -pRTv[ + pRv[T' 


P"v'U « o 


- _ - P" v k 

v* = V k - v k = IT 5 - - IT P,k 

P Pl°p 


V >"k + V"°ik ~ P{Tij,k + T ik ,j + Tjk.i) 


where 

k = ^r,, , C t = 0.11 [1 - exp(— y + /25)] , a f = 0.5 

Li 

In addition, we will use the model for II, y recently developed by Speziale, Sarkar and 
Gatski (1989), 

n,y * ?[(C lC + C 5 P)b {j - C 2 c(b ik b k] - |/J*y) 

- Csk{b ik Sj k + bjtSik — -b mn S mn 6ij) 

- C 4 k{b ik W jk + b jk W ik ) - 1(1 - C* II l ! 2 )kSij] 

o 
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where we have used the following notation 


P = > b a = u~ \ 6ij ’ 11 = bijbij 

Sij = 2 ■+■ Vj,i) > Wi j = — (Vfj — Vj i) 

The model constants for the term II<y are as follows: 

Ci = 3.4 , C 2 = 4.2 , C z = 1.25 

C 4 = 0.40 , C$ = 1.80 , C* = 1.62 

Finally, the model for the dissipation rate e,y is given by 

2_ x 

U, = 3 ptou 

Turbulent dissipation rate transport equation: 

The modelled transport equation for the dissipation rate e is as follows 

1 2 
dt(pe) + {pv k e),k = ^[-C e iPTij{vij - -v k ,kf>ij) - C f0 t> Jp,,] - C^p— 

— -Jfevk,k + {C t —Tkit,i),k + ( mc . i ) ,j ( 5 ) 

o € 

The model parameters in (5) are the following: 

C t0 = 1.0 , C tl = 1.44 

~pk? 

C t 2 = 1.83[l - exp(-iEy)] where R T = 

/X6 

C t = 0.15[1 - exp(— y + /25)] 
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Incompressible (Isothermal) limit: 


The Reynolds stress transport and dissipation rate equations reduce to 


d t {Tij) + (vkTij)'k = 
+ 


d t (e) + (wjfec),* = 


-TikVj t k - rjtVijk ~ II ij - -eSij 
2 k? 

g + Tik t j + 7y*,i)].k + [ 1/ { T ij,k + Tik,i + r jk,t)],k 

6 h 

—C e i—TijVij — C £2 — + (C e —r k ie i i) i k + 


where 


n?y = (Cie + CsP^j - C 2 e{b ilc b kj - hlS^) 

2 

^Z^{pik^jk "I - bjk^ik -bmnSmnfiij) 

- C 4 k{b ilc W jk + b jk W ik ) - 1(1 - C* II 1 / 2 )kS i j 

5 

and the coefficients are as defined before with y + = u T y/u. 


7 


NASA 

Naionai Aeronautics ana 
Scace Aamimstraton 

1. Report No. 

NASA CR- 18 1875 
ICASE Interim Report 8 
4. Title and Subtitle 


Report Documentation Page 


2. Government Accession No. 


3. Recipient's Catalog No. 


5. Report Date 


A Preliminary Compressible Second-Order Closure Model July 1989 

For High Speed Flows 6. Performing Organization Code 


7. Author(s) 


8. Performing Organization Report No. 


Charles G. Speziale 
Sutanu Sarkar 


Interim Report 8 

10. Work Unit No. 


9. Performing Organization Name and Address 
Institute for Computer, Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23665-5225 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23663-5225 

15. Supplementary Notes 

Langley Technical Monitor: 

Richard W. Barnwell 


505-90-21-01 

11. Contract or Grant No. 

NAS1-186Q5 

13. Type of Report and Period Covered 

Contractor Report 

14. Sponsoring Agency Code 


Final Report 
16. Abstract 

A preliminary version of a compressible second-order closure model that 
has been developed in connecton with the National Aero-Space Plane Project is 
presented. The model requires the solution of transport equations for the 
Favre-averaged Reynolds stress tensor and dissipaton rate. Gradient transport 
hypotheses are used for the Reynolds heat flux, mass flux, and turbulent diffu- 
sion terms. Some brief remarks are made about the direction of future research 
to generalize the model. 


17. Key Words (Suggested by Author(s)) 

Compressible Turbulence 
Second-Order Closure Model 
National Aero-Space Plane 


18. Distribution Statement 

34 - Fluid Mechanics & Heat Transfer 
Unclassified - Unlimited 

19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of pages 

22. Price 

Unclassified 

Unclassified 


8 

A0 2 


NASA FORM 1626 OCT 86 


NASA-Langley, 1989 



